NASA CONTRACTOR 
REPORT 


>© 


TT 





A CURVATURE-CORRECTED KIRCHOFF 
FORMULATION FOR RADAR SEA-RETURN 
FROM THE NEAR VERTICAL 


by Frederick C. Jackson 


Prepared by 

NEW YORK UNIVERSITY 
for Langley Research Center 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION • WASHINGTON, D. C. • APRIL 


2406 



1974 




1. Report No. 2. Government Accession No. 

NASA CR-2406 


4. Title and Subtitle 

A CURVATURE-CORRECTED KIRCHOFF FORMULATION FOR RADAR SEA-RETURN 
FROM THE NEAR VERTICAL 


3. Recipient's Catalog No. 


5. Report Date 
APRIL 1974 


6. Performing Organization Code 


7. Author(s) 

Frederick C. Jackson 


8. Performing Organization Report No. 

TR-72-1 


10. Work Unit No. 


9. Performing Organization Name and Address ^ ; , . 

New York University * * ' * f f ? Q. j 1 i .1 i 

School of Engineering and Science Contract or Grant No. 

Department of Meteorology and Oceanography NAS1-10090 


12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, DC 20546 


NAS1-10090 


13. Type of Report and Period Covered 

Contractor Report 


14. Sponsoring Agency Code 



16. Abstract 

A new theoretical treatment of the problem of electromagnetic -wave scattering from a randomly- 
rough surface is given. A high frequency correction to the Kirchoff approximation is derived 
from a field integral equation for a perfectly conducting surface. The correction, which accounts 
for the effect of local surface curvature, is seen to he identical with an asymptotic form found 
by Fock (1945 ) for diffraction by a paraboloid. The corrected boundary values are substituted 
into the far-field Stratton-Chu integral, and average backscattered powers are computed assuming 

-4 

the scattering surface is a homogeneous Gaussian process. Preliminary calculations for a K 
ocean wave spectrum indicate a reasonable modelling of polarization effects near the vertical, 

0 < 45°. Correspondence with the results of small perturbation theory is shown. 


17. Key Words (Suggested by Author(s)) 


Radar sea-return 


Curvature correction 


18. Distribution Statement 


Unclassified - Unlimited 


19. Security Oassif. (of this report^ 

Unclassified 


20. Security Classif. (of this page) 

Unclassified 


21. No. of Pages J 22. Price* 


For sale by the National Technical Information Service, Springfield, Virginia 22151 






















TABLE OF CONTENTS 


Page 


Summary v 

I. Introduction . 1 

1 . 1 Origins 1 

1.2 Rough-surface scattering theory and radar 

sea- return . 1 

1.3 Outline 5 

II. The Randomized Kirchoff Method According to 

Beckmann 6 

III. On the Nature of the Kirchoff Solution. Radar-Sea 

Return 15 

3.1 The stationary phase approximation 15 

3.2 The "small perturbation" approximation 19 

3.3 Radar sea-return. The idea of a composite 

surface . 21 


IV. Results of Small Perturbation Theory 27 

V. A High Frequency Correction to the Kirchoff 

Approximation * 30 

VI. Calculation of the Scattered Power 40 


VII. 


Properties of the Solution 

7.1 Comparison with first-order Rayle igh-R ice 
theory * 


7.2 A sample calculation: Copolarized returns 
from an isotropic RT^ spectral-law surface 
and a one-dimensional spectral-law surface . 


51 


51 


52 


VIII. Discussion and Conclusions to be Drawn, 


60 


References 


65 


Appendix 


69 


ili 



LIST OF FIGURES 


Page 

Fig. No. 

1 The scattering situation 7 

2 Defining geometry and notation „ . 7 

3 Two-dimensional Kirchoff solution for a K"^ spectral- 

law surface 24 

4 Stationary phase and small perturbation approximations 

to the Kirchoff integral 26 

5 The ’'composite surface" solution as sum of 

stationary phase and small perturbation cross-sections ... 26 

6 First-order Rayleigh-R ice solution 29 

7 Three-dimensional corrected Kirchoff solution for co- 
polarized returns from an isotropic K“ spectral-law 

surface 56 

8 Two-dimensional corrected Kirchoff solution for co- 
polarized returns from a s pectral-law surface 57 

9 Ocean backscatter data showing cross-over of copolarized 

returns .. 58 


xv 



Summary 


A new theoretical treatment of the problem of electromagnetic 
wave scattering from a randomly- rough surface is given. A high fre- 
quency correction to the Kirchoff approximation is derived from a 
field integral equation for a perfectly conducting surface. The cor- 
rection is of the for m = S where is the Kirchoff value 

of the current density, and S is a linear function of the second de- 
rivatives of surface height. The correction is seen to be identical'to 
an asymptotic form found by Fock (1945) in his investigation of dif- 
fraction by a convex paraboloid. 

The corrected current density is substituted into the far- 
field Stratton-Chu integral, and average backs cattered powers for 
the four linear polarization combinations are computed on the as- 
sumption that the scattering surface is describable as a homogeneous 
Gaussian random process. 

The strength of the solution is that local diffraction effects 
(arising from surface curvature) are properly correlated with surface 
height and slope without requiring their smallness. 

Application to radar backscatter and natural microwave 
emission from the sea is discussed. It is concluded that this "cor- 
rected Kirchoff" formulation offers a superior predictive capacity for 
co-polarized and depolarized returns from the near vertical, 0 < 45°. 
Extension to the bistatic case is recommended for application to the 
natural emission problem. The evolution of a new "composite" model 
of backscatter combining this solution with small perturbation results 
suggests itself. 


v 



I. Introduction 


1 . 1 Origins 

A satellite -borne combination radar-radiometer has been 
proposed as a remote sensing device for monitoring wind field/wave 
field conditions over the world ocean (see, for a recent version of 
this proposal, Moore and Pierson, 1971). This proposal has stimu- 
lated much of the recent theoretical and experimental work on radar 
backscatter and natural microwave emission from a w ind- roughened 
sea. In particular, the research reported here was motivated by a 
need for a broader theoretical basis for understanding and predict- 
ing microwave scattering by the sea. A good part of this research 
has been previously documented in a New York University technical 
report (Jackson, 1971). 

1.2 Rough surface scattering theory a nd radar sea-return 

The general problem of electromagnetic-wave scattering 
from a plane rough surface has been approached in a variety of ways. 
Some authors have dealt with exact solutions for scattering from 
certain simple (idealized) surfaces such as a surface composed of 
periodic rectangular corrugations (Deryugin, I960). Exact solutions 
for arbitrary surfaces can be obtained by the numerical solution of 
an integral equation for the surface field (Lenz, 1971). Except for 
one -dimens ional surfaces possessing a small amount of structure 
(e. g. , a small number of 1 ’hills" and "valleys") the computational 
time is prohibitive. Many analytical methods have been developed 
for the class of surfaces which can be called slightly rough, defined 
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as having small slopes and heights small compared to the wavelength. 
Exact solutions for deterministic slightly rough surfaces can be ob- 
tained by Rayleigh's method, or Meecham's (1956) variational tech- 
nique. Rayleigh's method has been randomized by Rice (1951), and 

this provides one of the most powerful methods for handling scatter - 

/ 

ing from randomly rough surfaces. Twersky (1957) has used a Ray- 
leigh image method for computing exactly the scattering by a random 
array of "bosses" or protuberances on a plane. Exact methods-- 
while impractical for application to randomly rough surfaces with 
a high degree of "structure"- -can be very useful in the testing of 
approximate methods developed for general classes of randomly 
rough surfaces. 

Some authors have taken conceptual approaches which 
are at the same time simple and instructive --for example, Katzin’s 
(1957) slope-facet model and Long's (1965) dipole model of sea re- 
turn. Katzin's model is interesting, for it entertains an important 
property of radar sea-return, that near vertical incidence (radar 
pointing downward) the backscatter mechanism is dominantly 
specular reflection, while backscatter from large angles of incidence 
is controlled by diffraction processes. Katzin's model has been ex- 
tended by Rouse (1970). 

Of approximate methods which have found practical appli- 
cation to scattering from continuously distributed random rough 
surfaces, there are essentially three: 

a) Geometrical optics 

b) Physical optics (Kirchoff theory) 
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c) Small perturbation (Rayle igh-R ice theory). 

Geometrical optics (or ray optics) formulations, because of 
their great simplification of the electromagnetic problem, are cap- 
able of handling such phenomena as shadowing and multiple scatter- 
ing (Lynch and Wagner, 1970). In considering the generally mild 
slopes of ocean surface waves, a ray-optical type of multiple scatter- 
ing is not likely to be a significant part of the backscatter mechanism. 
In the general bistatic case, however, the process of multiple 
reflections may be important. Shadowing effects can similarly be 
ignored for angles removed from grazing incidence. Shadowing 
is to a first approximation a simple enough process that it can be 
included in a physical optics formulation (Sancer, 1969). 

Most of the work in the last decade has been based upon 
either (or both) physical optics (randomized Kirchoff method) or 
methods of small perturbation. The randomized Kirchoff method 
developed by Beckmann (Beckmann and Spizzichino, 1963, Chs. 3 
and 5) and others uses the physical optics (Helmholtz) integral 
with the so-called Kirchoff or tangent-plane approximation to the 
boundary values of the field. The Kirchoff method is good for softly 
undulating surfaces having everywhere a local radius of curvature 
large compared to the electromagnetic wavelength. An advantage 
of the Kirchoff method (over small perturbation methods) is that 
surface height variations do not necessarily have to be small com- 
pared to the radar wavelength. A shortcoming of Kirchoff 
theory is that with its tangent-plane approximation it cannot account 
for polarization effects. Apart from its inability to account for 



polarization effects, Kirchoff theory has suffered because its in- 
herent strength was not exploited. Often, what amounts to a station- 
ary phase approximation to the Kirchoff integral is made. The 
stationary phase approximation is equivalent to geometrical optics; 
so, the ability to account for diffraction effects is lost. Chia (1968), 
in applying Kirchoff theory to radar sea-return, appears to be the 
first to have avoided this approximation by using a realistic wave- 
height covariance function in the Kirchoff integral. 

Small perturbation theory has gained increasing favor in 
the last few years among scientists working on radar sea-return. 

This is primarily because of its ability to account for polarization 
effects but also because of its simplicity. The explicit dependence 
on the wave-height spectrum pointed the way to using a realistic 
representation of the rough sea surface (Valenzuela, 1968; Wright, 
1968). Small perturbation methods have been developed by several 
authors (Bass, 1961; Wright, 1966), notably by R ice (1951). Rice^ 
randomized Rayleigh method seems to be the superior, for it is capable 
of iteration to higher order in height and slope, and is capable of pre- 
dicting depolarization in the plane of incidence (Valenzuela, 196 7, 

1968). The Rayleigh-Rice method- -although having some com- 
monality with the Kirchoff method--is a fundamentally different 
approach to the scattering problem; whereas Kirchoff theory is a 
high frequency treatment, Rice's theory is a low frequency method. 

In the high frequency approximation, the field in the vicinity of a 
surface point is dependent only on the local geometry of the surface, 
i. e. , height and slope. In the Rayleigh approach, the field at a point 
is related to integral properties (rather than local properties) of the 
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surface. This local versus integral (or modal) duality in electro- 
magnetic theory is discussed by Felsen (1964) in his review of high 
frequency diffraction. 

1 . 3 Outline 

The theoretical approach taken in this work is a high- 
frequency one , and is essentially an extension of the randomized 
Kirchoff method. A curvature correction to the Kirchoff approxi- 
mation is derived from a field integral equation, and the corrected 
boundary values are used in the Stratton-Chu far-field integral 
(vector form of the Helmholtz integral). 

In the following section, Beckmann's development of scalar 
Kirchoff theory is recapitulated in order to place the theory of 
Sections V and VI in its proper perspective. Section III contains 

a discussion of the Kirchoff solution, and presents the results 

- 3 

of Kirchoff theory applied to a scattering surface with a K spectral 
law. Section IV presents the results of small perturbation theory. 
Because of the complexity of the scattering integrals arrived at 
and their requirement of a detailed knowledge of the wave-height 
covariance function, no thorough computation and comparison with 
sea-return data is made. However, some sample calculations are 
given, the nature of the solution is discussed and estimates of its 
strength are made. 



II. The Randomized Kirchoff Method According to B eckmann^ 

Consider the scattering situation depicted in Fig . 1. 

A monochromatic source (radar) illuminates a portion S o of a rough 

conducting surface. Beckmann, deals with the scalar electric field 

E scattered by the surface. The surface field vanishes outside of 

S . and the scattered field satisfies the radiation condition at 
o 

infinity. The electric field scattered toward an observation point 
P is given by the Helmholtz integral involving the (unknown) field 
boundary values: 


E(P) 




E 


9G 

9n 


r £E\ 

G an ) ds 


( 2 . 1 ) 


where E is the total electric field or. the surface; 0/3n stands 
for the normal derivative directed outwards (upwards) from the 
conducting volume; and G is the Green's function for (three- 
dimensional) free space, 




* ikr 


where k is the microwave propagation constant and r is the 
distance from source points on the surface to the observation 
point. If P is in the far-field (Fraunhofer zone) r can be 


*In Beckmann and Spizzichino (1963), Chapters 3 and 5. 
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Conducting Volume 


* Fig. 1. The scattering situation. 



Fig. 2. Defining geometry and notation. 
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approximated by' ,: 

r ~ R - R • x (2.3) 

where R = RR is the position vector of P and x is the coordinate 
vector of surface source points. With this approximation, we can write 

- ikR LkR • x 

G=| ' 6 (2.4a) 

and 

|- = ikR-nG (2.4b) 


"'Here, Beckmann fails to give any criterion for deciding how far 
removed the observation point must be to lie within the Fraunhofer 
diffraction zone. A criterion can be established, however, without 
too much difficulty. The phase error incurred in making the far- 
field approximation is (see accompanying figure): 



A<p = ~ ( XP - QP ) 


2tt 1 2 

X • 2 a • 


With a ~ x cos 9/R, and 
setting QP ~ OP = R 




2 

TTX COS 0 

XR 


The first Fresnel zone is at = tt, so that \f\R sec 0, 

For example, with a 10 cm radar at an (aircraft) altitude of 
1 km, 10 m = 100\; or, at a (spacecraft) altitude of 1000 km, 
x ^ 100 m = 1000\. For deep fade conditions (rms surface height 
much larger than radar wavelength), the surface source field be- 
comes incoherent well within one -hundred radar wavelengths. 
Only under laboratory conditions and/or the condition of small 
roughness amplitude need Fresnel zone effects be considered 
(see, for example, Wright and Keller, 1971). 
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l 

where n Is the unit outward surface normal. 

The Incident field E L is taken to be a plane wave ' of unit 
amplitude. With the time -harmonic dependence e La ^ suppressed, E l 
is written as the phasor, 

i - ik • x 
E = e • — 1 — 

where k^ is the propagation vector of the incident wave. 

The Kirchoff approximation to the field boundary values E 
and 3E/9n consists in assuming that the field in the vicinity of (an 
"epsilon" neighborhood) and at a point on the surface is nearly equal 
to the field which would exist on an infinite tangent plane at the point. 
This is a type of "high frequency" approximation. For a high enough 
microwave frequency (wave number), the surface curvature "appears 
mild to the radiation and the surface can be considered to be locally 
flat. There is then a perfect reflection of the incident wave in 
accordance with the geometrical optics Law of Reflection; the ampli- 
tude and phase of the reflected wave are given by Fresnel’s formulas. 
Then, writing the field in the vicinity of the surface as an incident 
wave plus a reflected wave, we have for the total field and its normal 
derivative evaluated on the surface: 

E = (1 + ftjE 1 

3E ~ . _ i * ** 

9 ^ = - ikj -n(l - R)E . 

* The justification for ignoring the sphericity of an incident wave 
radiated by a point source (radar) follows the same arguments 
given for the far -zone approximation on the scattered field (see 
previous footnote). 

**In writing (2.6) with the incident wave E 1 given by eq. (2.5), it is 
tacitly assumed that there is no multiple scattering nor shadowing 


(2.5) 


(2.6a) 

(2.6b) 
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The Fresnel reflection coefficient ft is in general a matrix. Treat- 
ing ft as a scalar coefficient might seem to be of dubious validity, 
but it will serve us in a formal development. 

Brekhovskikh '' has shown that the Kirchoff approximation 
is valid if 

47rr c co s 0 1 >> X (2.7a) 

where r^ is the radius of curvature, 0 1 is the u local M angle of 
incidence (the angle included between the local normal n and the 
incident ray), and X is the radar wavelength. Wait and Conda v * 
have given the criterion, 

■rrr c cos^G 1 >> X . (2.7b) 

These inequalities' should not be interpreted in a strict sense, since 
large third derivatives can exist even if the second derivatives are 
small. We should really interpret r^ as a root-mean-square value 
for the surface, giving an indication of the degree of "smoothness 11 . 

Substituting eqs. (2. 4-2. 6) into the Helmholtz integral (2.1), 
we get for the far-zone scattered field: 

-ikR i(kR-k ) • x 

E=—^~ / [ikR.n(l+ft)+ ikj-nfl- R)]e ■ 1 dS (2.8) 

S o 

Departing from Beckmann's' development, we specialize 
the problem to the monostatic case (backscatter). The observation 

'See Beckmann and Spizzichino (1963), p. 29. 

** Ibid . , p. 29. 



point then coincides with the source, and the incident propagation 
vector can be written as kj = -kR . Eq. (2.8) becomes: 


E = 


ike ~ i2kR . x 

/ 2& R • n e — dS 

S 

o 


4ttR 


( 2 . 


We adopt the following notation (Fig. 2): Let the surface be 

described in the (x, y, z) Cartesian system by 


z = f(x, y) . 

The coordinate vector of surface points is then 

x = (x, y, f) 


The normal vector is given by 


n = (-f x> -f , 1 ) cosw , 

cos w = (1 + f 2 + f Z ) 2 
' x y 

where f and f are partial derivatives. The surface area element 
x y 

is given by 


dS = sec co dx dy . 

The plane of incidence is formed by R and the z-axis. The x and y 
axes are oriented so~that the x-axis lies in the plane of incidence. The 
angle of incidence 0 is measured positive where the incident ray comes 
from the left (negative x-direction). The unit vector R is given by 

R = (-a, 0, y) 


whe re 


a = ' s in 0 
y = cos 0 


Equation (2.9) becomes 
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E= . //2R- (v + af ) •e' i2k(aX '^ f, dx dy (2.10) 

A 

o 

where A q is the (horizontal) illuminated area. Now assume for 
the moment that the surface is perfectly conducting. Then, ft = ±1 
(the sign depending on polarization) and ft can be treated as a con- 
stant and removed from under the integral/' With ft = ±1, (2.10) 

is integrated by parts to yield 

’“// e " ^k( ax ~Vf)dx v dy + n edge terms M . (2-11) 

A o 

2 

For large areas k A q >> 1, the edge terms are negligible. This 

is not obvious from the form of (2.11); but, it is physically reasonable 

that-provided the surface is moderately rough and grazing incidence 

is avoided --edge effects are unimportant in the scattering problem. 

2 

The average backs cattered power is proportional to ( |E| ) where the 

brackets denote expectation, or ensemble average. Neglecting the 

2 

edge terms, we form |E| = EE* as a two-fold integral over A q , and take 

the average over all realizations of the surface: 

< |E|2)= (j2^) 2 ////< ei2kY(f, " f) > e “ i2kQ(X '' X) dx ,d y , dxdy (2.12) 

The expectation 

r 2k^y) = (2.13) 

is the two-dimensional characteristic function of the random vector 
(f ! , f) evaluated at (2ky, - 2ky). If f is a stationary (homogeneous) 

*This is not immediately evident. It turns out, however, that the 
vector formulation of the Kirchoff integral for perfect conductivity 
yields precisely (2.10) with the effect of polarization properly accounted 
for by the ± sign in front of the integral. (Cf. eqs. 6.2a, d) 
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Gaussian random process of zero mean, then 

„ = e 4k\ 2 B(i) 

where 

B(t) = -R(0) + R (|) 

and R is the surface height covariance function 


(2.14) 


(2.15) 


R(|)=<f«f> (2.16) 

and ^ is the "lag n or ’’separation" vector, 

1 = r)) = (x 1 - x, y f - y) . (2. 1 7) 

2 

In terms of the surface height variance, o* = R(0), and the auto- 

2 

correlation coefficient, p = R(£)/<r , B can be written as 

B(|) - - <r 2 + R(|) = - cr 2 (l - p) . (2.18) 


If is large compared to the scale of roughness corre- 
lation lengths") the double-area integral (2.12) is nearly equal to 


, . t-, . 2 / k \ 2 °f 7 4k 2 Y 2 B -i.2ka£ 

<' E ' >=2^/0/ / e Y e 

‘ / - OO - OO 


d| dr\ 


(2.19) 


The infinite limits on the r\ integration are not to be taken in a strict 


sense. For large r\ the exponential approaches the constant asymp- 

2 2 2 

totic value exp{-4k y cr } . The integral over r\ then increases 


linearly with rj , and is independent of £ . If the infinite r| limits 
are interpreted in the sense of the limit 
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llm 
Y — oo 


7 T 4k v 2 B -i2ka£ JC , 

/ / e * e dT] 

-oo -Y 


then it is seen that (except for vertical.incidence, a = 0) the 

phasor nullifies the constant contribution from large r\ and the integral 

exists. In practice we generally deal with very rough surfaces for 
2 2 2 

which 4k y o' >> 1 and the exponential for large lag is extremely 
small. A numerical integration over n would then be stopped, at some 
prescribed value of r| for which the exponential is very small. 

Since the covariance function R(4) is symmetric about 
the origin, (2.19) simplifies to 

2 oo oo 2 2 

< aJ cos 2ka| / e4k Y Bdr ' • < 2 * 2o > 

L 71 V O - oo 


The backscattered power is conventionally given in terms of the 
normalized isotropic radar cross section, cr°, defined by 

2 

, <|E| 2 > . (2.21) 

o 

Then, (2.20) becomes 

7 ,2 oo oo . 2 2 

o-° = ~ • / /e v a cos 2ka| d£ dt} . (2.22) 

y O - oo 
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III. On the Nature of the Kirchoff Solution. Radar Sea-Return 

The Kirchoff Integral (2.22) contains information on the two 
backscatter regimes: the ’'specular" (ray-optic) regime which pre- 

ponderates for small angles of incidence, and the "small perturbation" 
(diffraction) regime which accounts for the backscatter from large 
angles of incidence. 

3. 1 The stationary phase approximation 

The high frequency limit of physical optics yields geometrical 

optics : symbolically, 

lim (physical optics) = geometrical optics 
k—**oo 

The high frequency limit is equivalent to a case where (1) the surface is 
very smooth, and (2) the phase modulation is very deep. The "deep fade" 
condition ( Hagfor s f 1 964 ) means we must have (for backscatter): 

4k 2 v 2 (T 2 >> 1 . (3.1) 

In the high frequency case characterized by (3.1), the exponential 
2 2 

exp{4k y B} becomes negligible outside of a small region about the 
origin. The meaning of the smoothness condition is that within this 
"neighborhood" of the origin, the covariance function can be repre- 
sented by a Taylor series truncated at second degree. Thus, we can 
write B(_|) as the paraboloid: 

m2 2 

B = R (0) 1 + R (0)£r) + R (0) \ 

XX- 2 xy — 1 yy'— 7 2 

Considering (for simplicity) the special case where £ an< ^ are the 
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principal axes of the ellipse B = constant, B can be written as 

R 1 2,2 i 2 2 . 

B = " 2 % ' 1 a y ’l < 3 ‘ 2 > 

2 2 

where cr = -R (0) and cr = -R (0) are the slope variances in 
x xx'~ y yy — r 

the x and y directions. The rms slope cr is invariant with re- 
. s 

spect to coordinate rotation. 


cr = 


s 




Putting (3.2) into (2.22) and integrating, we find (letting o^ 0 
represent the high frequency limit): 



2<r cr v 
x y Y 


1 tan^9 

2 , 2 
x 


(3.3) 


(3.4) 


The cross-section (3.4) is independent of k, consistent with the fact 
that it is the high frequency limit- -geometrical optics. Another way 
of seeing that (3.4) is a geometrical optics limit is to note the cr 0 
is simply proportional to the probability of a surface facet satisfy- 
ing the specular (ray-optic) condition for backscatter: = tan0. 

For an isotropic surface, cr x = cr , and (3.4) becomes 

t 2 . 
tan 0 


1 -s 2 

I""? e 


(3.5) 


<r v 
s Y 


For reference, we should like to know the cross-section for a two- 
dimensional scattering situation. One has to go back to the Helm- 
holtz integral for two dimensions: 
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E= ?i.)( E lf - G lf ) dS • 

So 

where G is the zeroth-order Hankel function of the second kind: 

G = H^ 2) (kr) . 

Making the far-zone approximation and following a procedure identi- 
cal to the three-dimensional case, but defining the cross-section 
for cylindrical waves as 


<r° = <|E| 2 > 


(3.6) 


where X is the linear extent of illumination, we get: 


o* 


oo 2 2 

2k r 4k y B 

2 J e 


cos 2ka£ d£ 


(3.7) 


The high-frequency limit of (3.7) is: 


2 

i. t an 9 

_ o _ n/tt 7 2 ^ ^x 2 

°1 ~ 3 e 

V 


(3.8) 


Taking the high frequency limit of the cross-section (2.22) 

is seen to be equivalent to making a ’’stationary phase 1 * approximation ' 

to the field integral (2.10). The phase <j> = 2k(ax- yf) is "stationary 11 

for the specular points having slope f = a/y = tan 0 and f =0. 

x y 

The stationary phase approximation provides a (mathematical) 
rationale for accounting for finite conductivity effects in the specular 


"'For example, see Lamb, Sir H. : Hydrodynamics. Dover Publications, 
New York, 1963, p. 395. 
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(near vertical) regime. Referring again to eq. (2.10), we see that 
ft is brought out of the integral and evaluated at the local angle of 

t 

incidence 0 = 0 for stationary or specular points. This is just 

the mathematical version of what we should expect on physical 

grounds: that where specular reflection is the dominant mode of 

scattering, the backscattered field is composed primarily of incident 

waves reflected at locally vertical incidence. Thus, the effect of finite 

conductivity it to a first approximation (a very good first approximation 

for radar return from the sea) simply to reduce the backscattered 

2 

signal by a factor of |ft(0)| . And this is true regardless of polari- 
zation, since at (locally) vertical incidence the distinction between 
horizontal and vertical polarization' 1 ' disappears. 

Stogryn (1967a) came to this conclusion using a vector 
Kirchoff formulation and making the "stationary phase 11 approximation. 
Kaufman (1 971 )- -apparently unaware of Stogryn ! s conclus ions- -used 
a vector Kirchoff formulation for finite conductivity, and computed 
cross sections for different polarizations. Unfortunately, Kaufman 
made the stationary phase approximation, so that no information was 
gained on the effects of finite conductivity (as they manifest in a 
Kirchoff formulation) in the diffraction regime at larger incidence 
angles. For reference, let us write down our conclusion symbolically 
as 


*These terms have not yet been defined. Vertical polarization (V) 
means the E -vector is in the plane of incidence (mean plane or 
local plane); horizontal polarization (H) means the E -vector is per- 
pendicular to the plane of incidence, that is, lies in the plane of 
the surface (mean or local). 



(3. 9) 


*{ = 1^(0) | 2 (Tj- <00) 

where a*°(oo) stands for the cross-section of a perfectly conducting 
surface, 

3.2 The n small perturbat ion M approximation 

Small perturbation theory has been remarkably successful 

in predicting sea-return from large incidence angles (e. g. , see 

Guinard et al. , 1971). The reason is that toward larger angles of 

mcidence a larger portion of the surface becomes effectively small 

amplitude by the Rayleigh criterion. The Rayleigh criterion is actually 

2 2 2 

contained in the Kirchoff formulation as the exponent 4k y cr ( cf . 
Beckmann and Spizzichino, 1963, p. 93). A surface is smooth by the 
Rayleigh criterion if 

2k^o~ << 1 . (3.10) 


Although the sea-surface variance is generally several 

orders of magnitude too large to satisfy (3.10), what happens is that 

2 2 

at larger angles, the exponent 4k y B remains small over larger lags. 
Also, the horizontal wave number 2ka increases, causing relatively 
more energy to be contributed (to the Kirchoff integral) from the 
small lag neighborhood where the exponent is small. 

Assuming that (3.10) is satisfied, we can approximate the 
exponential by its first-order expansion, 


4k 2 Y 2 B 
e Y 


1 - 4k 2 Y 2 cr 2 + 4k 2 Y 2 R(|) 




(3. 11) 
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The Kirchoff Integral (2.22) then becomes, approximately ( naming 
the small perturbation cross-section o'^): 

4 oo oo 

o " 2 = — — — / / ^ (.i) cos 2ka£ d£ 

o - oo 

+ a specular term at 0 = 0 (a Dirac spike) . 

From the definition of the wave-height (roughness-amplitude) spectrum 
(half-spaced): 

S(j_) = JR(S) cos(£- t) d£ dr, (3.12) 

TT O - OO 

we see that 0*^° can bo expressed simply as 

<r 2 ° = 8tt k 4 S(2ka, 0) . (3. 13) 

The wave number = (2ka, 0) is referred to as the Bragg back- 
s catter condition (Wright, 1968). The physical meaning is understood 
as follows: Huygens wavelets are in phase in the backscatter di- 

rection when the source elements are located a distance 

L = y esc 0 

apart. And we have i - 2tt/L = 2ka. 

We should also like to know the two -dimensional form of cr^ . If 
the (one-sided) one-dimensional spectrum is defined as 

s(!) = §/ R(6) co S {| de , (3.14) 

o 

then applying (3. 11) to (3.7) we get: 

<7 2 ° = 4 k\s(2ka) . 


(3.15) 



3 . 3 Radar sea-return. The idea of a composite surface . 


Chia (1968) appears to be the first to have used a realistic 
wave height covariance function in a Kirchoff formulation. What Chia 
did, essentially, was to take the equilibrium range spectral law 
(Phillips, 1966) and cosine transform to find the covariance function. 
The equilibrium range spectral law is (in polar coordinates): 

S(K,cp) = AK"^F(cp) (3.1 

where A is a universal constant and F is a dimensionless spreading 
factor normalized so that 

Zn 

f F(cp) dcp = 1 . 
o 

The equilibrium range is supposed to exist in a fully aroused sea be- 
tween wave numbers near the spectral peak down to wave numbers 
approaching the capillary-wave regime. Assume that isotropic 
conditions prevail, and that the spectrum can be defined by 


A -4 

S(K) = ~ K 
t r 


|cp | ^ TT 


K ^K< 00 
o 


(3.1 


where K q is a low wave number near the spectral peak. 


K 0 ~ g/U 2 ; (3.1 

where g is the acceleration of gravity and U is the wind speed at 
a nominal anemometer height. 

For microwaves, the propagation constant k is several 
orders of magnitude larger than K q . Scattering is determined pri- 
marily by the wave structure in a wave number domain centered about 
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the Bragg backscatter condition, K = 2ka. This allows a simple ap- 
proximation to the covariance function corresponding to (3.17). For small 

lags K o r « 1 = j ? + r\ 

A °° 11/2 -3 

R(r) =— f / K cos(Kr coscp)dK dcp 
K o -” /2 

2 2 

~ (T + A ~ (-1 + y + In K r / 2). (3.19) 

o 

where y ~ 0.577. . . is Euler’s number. 

Since the high-frequency portion of the wave spectrum 
Is always nearly Isotropic, a one -dimens ional counterpart to (3.16) 
can only be a fictitious analog. However, if one imagines all the 
wave energy to be concentrated into one direction, then F = 6(cp), 
a Dirac spike, and the one-dimensional spectrum becomes 

S(K) = AK" 3 (3.20) 



The approximation to the covariance function for small lags, 

K q £ << 1 , is 

00 3 

R (|) = A f K cos K6 d| 

K o 

2 3 

- <r + A-f- (- -£-+ Y + In K q £ ) (3.21) 


It is possible that funnelling all the wave energy into one 
direction might produce an unrealistically large amount of scattering 
(i.e. , in the two-dimensional approximation to the three-dimensional 
scattering problem). Let us allow for a variation of the spectral 
constant A by imagining that only a certain band of directions 

*Luke, Y. L. , Integrals of Bessel Functions , McGraw-Hill, 1962, p. 48. 
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are funnelled into the one direction, that is, let -A —►A —A. If 

TT 

nothing else, this artifice will let us see the variation of <r° with a 
changing spectral constant/ 1 ' 

Figure 3 shows the two-dimensional Kirchoff solution 

(eq. (3.7)) using the covariance function given by eq. (3.21). Nominal 

-3 - 1 

values of A - 5 x 10 , U = 15 m sec and X. (radar) = tt cm were 

u sed in the calculation. The effect of varying the low wave number 
cutoff is small (on the order of a few decibels for wind speeds between 
5 and 20 m sec ^) and is not illustrated. The effect of increasing 
wind speed is simply to cause a greater ,, tilting ,, of the small wave 
structure (which is primarily responsible for the scattering) by the 
larger waves. The result is a "smearing 11 of the incoherent scattered 
power pattern. Also not shown is the frequency dependence of the re- 
turn. The frequency dependence can be expected to be small. This 
follows from the k- independence of both stationary phase and small 

3 

perturbation approximations. (Note that k S(2ka) is k- independent 
- 3 

for the K spectrum. ) 

A comparison of stationary phase and small perturbation 

approximations with the Kirchoff integral is made in Figure 4. A 

value of o- = 0.215 was chosen to match c r 0 with the Kirchoff <r° 
x 1 

at vertical incidence. 


'■'The spectral constant is wind-speed dependent. The observations of Ley- 
kin and Rosenberg (1970), for example, show that A increases 
with wind speed until for high wind speeds A levels off and 
approaches its asymptotic (equilibrium) value. 
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The coincidence of the small perturbation approximation 
(o^ 0 ) with the Kirchoff curve provides a validation of small perturb- 
ation theory for sea-return from larger incidence angles. 

The basic idea behind various "composite surface" models 
(e. g. , Semyonov, 1966; Fung and Chan, 1969; Krishen, 1971) is illus- 
trated in Fig. 5. The backscattered power from small perturbation 
and geometrical optics regimes is added incoherently. 


o 

c r 



( 3 . 22 ) 


where the specular (near vertical) portion of c r^° is suppressed. A 
remarkably close agreement with the full Kirchoff integral is achieved. 



in Decibel* 


Angle of Incidence 6 in Degrees 

Stationary phase and small perturbation approximations to 
the Kirchoff integral., The Kirchoff solution is the same as 
curve 'a 1 in Fig. 3. cr^ and are respectively given by 
eqs. (3„ 8) and (3. 15) 0 


Angle of Incidence 6 in Degrees 


The "composite surface" solution as a sum of stationary 
phase and small perturbation cross sections. 
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IV. Results of Small Perturbation Theory 

Methods of small perturbation apply to surfaces for which 

|f x l << 1 , | kf | << 1 . (4.1) 

The effect of surface roughness is then In the form of a small per- 
turbation on the primary field on a perfectly flat horizontal surface. 
Various small perturbation methods have been advanced; for ex- 
ample, by Rice (1951), Bass (1961), and Wright (1966). Rice used 
a randomized Rayleigh method. The field and surface height are ex- 
panded in a Fourier series. The surface height is expressed as 

f(x, y) = S P e i a (mx+ny) (4.2) 

J mn 

m,n 11111 

where the P mn are independent Gaussian random variables, and 
a = 2 t t/X is the fundamental wave number. First-order Rayleigh - 
Rice theory for perfect conductivity yields a far-zone cross section 
of the form 

<t° = 8trk 4 g(0) S(2ka, 0) (4.3) 

where the spectrum is defined by equation (3. 12) The angle function 

g(0) depends on the polarization. For vertical polarization (Ei-vector 

in the x-z plane), ? ? 

g vv (0) = (1 + s in 0) . (4.4a) 

For horizontal polarization (E 1 -vector parallel to y-axis); 

g HH ( 0 ) = cos 4 0 (4.4b) 

First-order Rayle igh-R ice theory gives no depolarization in the plane 
of incidence. To account for depolarization (cross -polar ization) , one 
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one must go to second order in the ordering parameters |kf| and 

| f x | • This has been done by Valenzuela (1967). Fig. 6 compares 

the Kirchoff and Rayleigh-R ice solutions for the surface described 
_3 

b y the K spectrum. 




I 1 I I I 

20 30 4 0 5 0 60 

Angle of Incidence 6 in Degrees 


-3 


A « 1.25 x 10 




30 


V. A High -F requency Correction to the Kirchoff Appro xim ation 


A simple correction to the Kirchoff approximation that 
accounts for the effect of surface curvature can be extracted from a 
field integral equation. For a perfectly conducting surface free from 
singular curves (cusps), the integral equation for the magnetic 
field is (e.g., Fock (1945)): 

x V r G dS' (5.1) 


J 

— s 


— s 


+ n X 


-kUi 


where J = n X H is the surface current density and H is the 
magnetic field on the surface; where: 

= n X ZH l is the Kirchoff value of the current density, 
n is the unit surface normal directed outward from the 
conducting volume, 

H l is the incident magnetic field, 

( 1 ) (prime) denotes source point coordinates x 1 as opposed 
to field point coordinates, x , 

G = exp(-ik| x - x 1 | )/|x - x l | is the Green's function for 
homogeneous space. 


The assumption of a first-order continuous surface is con- 
sistent with high frequency approximation we shall be making. 

That is, as with the Kirchoff approximation, we shall require gentle 
curvature (\/r <<: 0* Mittra has derived a set of field integral 

equations which are valid for surfaces possessing sharp edges. 
These equations may provide a basis for future research. The 


*Mittra, R. A. J. : Notes for a course given at the University of 
Illinois, Department of Electrical Engineering (unpublished). 



formulation (5.1) is ideally suited to near-planar or smooth geometries 
For, when the surface is locally nearly flat or nearly planar, the inte- 
gral contribution is small compared to the Kirchoff value n X 2H 1 and 
can be regarded as a f, small perturbation’ 1 on the Kirchoff value. This 

is seen from the fact that the vector J 1 X ^'G is oriented nearly 

— s 

pa rallel to the field point normal n, so that the cross-product 
n X J 1 X V'G is small. 

Fock (1945) has examined the conditions under which a 
perfect conductivity formulation is valid. An assumption of perfect 
conductivity is strictly justifiable only when the radius of curvature 
is large compared to the skin depth (depth of penetration of the field). 
Conduction currents can then be considered to be confined to a thin 
skin layer, and can then be represented by a surface current density. 
Thus, we must have 

6 << r 

c 

where 6 is the skin depth (normal to a flat interface), the depth at 
which the field has attenuated by a factor of e ^ . For all microwave 
frequencies up to and including X-band (3 cm), sea water is a good 
conductor in the sense that the skin depth is much smaller than the 
wavelength. Even at X-band, 6 is small: 6/\ ~ 1/2-rr. Thus, a 

perfect conductivity formulation of the integral equation is consistent 
with the high-frequency (small \/r ) approximation we are making. 

We develop the integral equation for the horizontal rough 
surface, z = f(x, y). We use the following notation: 

‘From Saxton and Lane's (1952) data. 



32 


x = (x, y, f), x* = (x 1 , y' , r 
p = x 1 - x = (£, q, L) 


n = f-f , - f , 1 ) cos gj 

: x y 


2 2 -- 

cos w - (1 + f + f ) *' 

x y 


dS = sec oj dx dy 
If we define the two-vector J as 

J = I I ~ | | sec gj 

then the three-vector equation (5. 1) can be written as the two- 
vector matrix equation, * 


mm mk 

J 


"x • J 

X 


— s 

J 


y-J 

_ y_ 


__ — s 


J = • J (0) + ^ ij; _M JJ d | d ri 


where 


-ikp 


^ = - (1 + ikp) — 


and M_ = M (x; £) is the matrix 


M = 




^ { y + z 


- n 


ax 


- e - 

5 ay 


, fi-er + t 

OT) X 


In this notation £ is to be considered a function of x and 

C = t, (x; £) . 

so that 


See Jackson (1971), pp. 4-8, for details. 


(5.2) 


(5.3) 


(5.4) 


(5.5) 
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§s- 


and 


d t 1 

= f . f 

3x x x 


(x — y) 


(x — y) 


The J component follows from the condition that J is tangential 
z — s 


to the surface: 


J = f J + f J 
z xx y y 


(5.6) 


We relax the high frequency (small curvature) condition of 
the Kirchoff approximation to allow for some degree of curvature. 

Since we are still dealing with a "smooth" surface, we can assume 

the surface height f has a Taylor series expansion about every point, x : 


f(x) = f + f u + f v + f - 5 - + f uv + f 3 - + • • * . 
x y xx 2 xy yy Z 

From now on it will be understood that f and the derivatives are to 

~ be evaluated at the local origin x^. We shall have u = (u, v) stand 

for the relative (horizontal) position vector of a field point. The 

relative position vector of a source point shall be given by 

1 - Xq = u + . If we form the difference £ = V - f we get 

,2 2 
C = f x S + f y 11 + f XX ( U £ + f-) + f X y( UT l + + £t)) + f yy (vr) +■£“) + ••• 


A little algebra will show that the matrix M has the expansion 

, 2 


M = M (x; %) » 


2 

-f 5- + f f- 

xx 2 yy 2 


•f r\ - f ip 
xy 1 xx ' 


-f i - f in 

xy yy* 


1 2 2 
f — f H- , 

XX 2 yy 2 J 


(5.7) 
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Thus, to a first approximation, M is independent of the local field 
point coordinates u and can be expressed in terms of the separation 
vector ^ alone. Call this approximation M^. 

A first approximation to \\) ls obtained by letting p cor- 
respond to the distance on the tangent plane, 

1 

P = P l = (aV + b 2 ^ 2 ) 2 (5.8) 

whe re 

2 2 2 Z 

a = 1 + i » b = 1 + f . . (5.9) 

x y ' 


If we let 


-ik Pl 

4, = = - (i + ik Pl )^-^ 

Pi 


(5.10) 


the leading error term is proportional to the curvature. 

If the approximation ^ M - is made the leading 

error term is proportional to third derivatives of surface height and 

the product of two second derivatives. What we are going to assume 

is that the bulk of the integral / dt] is formed in the 

neighborhood of x where the error terms are small. This as- 
® — o 

sumption is difficult to justify. Cullen (1958) has examined the high 
frequency behavior of the integral equation (5.1) for a convex body. 
There is an apparent contradiction in the relative importance of the 
source distribution near and removed from the field point. We will 
see this when we compare the radiative (far-source) contribution 
to the inductive (near - source) contribution to the integral (see below, 
p. 39). We should expect the mild curvature restriction on this 
approximation to be very similar to the criteria (inequalities (Z. 7)) 
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for the Kirchoff approximation, except that this approximation should 
2 2 

be good to 0(\ /r^ ) rather than 0(\/r c ) . 

In any event, assuming the integral 

equation (5.3) becomes 

J(u) = J (o) (u) + ~.j /4 j(1) (p 1 ) M (2) (f) dt! (5.11) 


Equation (5.11) can be solved exactly by Fourier transformation techni- 
ques. Fourier transformation is now a common method for solving 
two-dimensional (plane) diffraction problems, where integral equa- 
tions of this type occur (Bouwkamp, 1954). But, remember, unlike 
a true two-dimensional equation, eq. (5.11) is only approximately 
correct. The older method of iterated kernels is a more appropriate 
method of solution. 

The integral operation takes an 0(1) quantity into an 

2 2 

0(X./r ) quantity; and 0(\/r ) quantity into an 0(\ /r ) quantity, 
c c c 

and so on. Since J differs from by an 0(\/r c ) quantity, we 

can set 


J = J<°> + J (1) 

where = 0(1) (for a unit magnetic field) and = 0(x/r c ) . 

Equation (5.11) then yields for 


(5.12) 


J (1) (u) = — //q J (1) (p 1 )M (2) (i) J (0) (u + l)d| dq + 0(\ 2 /r c 2 ) . 

We have really lost nothing here since equation (5.11) was accurate 
only to Q(\/r^) to start with. Note that the above equation is 
most accurate at the local origin u. And since we no longer need 
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the convolution properties of the original equation because is 

a known function, all we need calculate is 





d| dr\ 


For a plane-wave incident field with unit magnetic vector, 
has the form 


J (o) = 2Ae" ik(aX -V f) 


where for vertical polarization (E l -vector in the x-z plane): 



and for horizontal polarization (E L -vector parallel to the y-axis): 

- af 

Y 

v + af 

Y x -* 



It is entirely consistent with the development of equation (5.13) to 

expand in terms of £ about Xq neglecting terms of 0(\/r ). 

We can do this because when multiplied by the kernel, 0(\/r ) 

2 2 

terms become 0(\ /r ) terms. We then have approximately 
J (o) = 

where A_ is equal to (5.15) with f^ and f evaluated at the local 
origin, x- (The sub-zero notation is abandoned.) And is 
the wave number, 

1 = k(a - Y f x > -\ f y ) 

Putting in the integral (5.13) we get 

jt 1 ) = § J (o) 


(5.13) 


(5.14) 


(5.15a) 


(5.15b) 


(5.16) 


(5.17) 


(5.18) 
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where S_ is the Fourier integral. 




(5.19) 


The infinite limits have been applied just to make the integral definite. 

The evaluation of this integral is given in the Appendix. 

The elements S.. of S are found to be: 

U 

,1- m 2 )% - (1 - ( 2 ) V 

b^J 


11 2kabv 3 


12 


’21 


f 


- i 

(i 



, i 3 

kabv 

- 

a 



f 

f n 

- i | 

■I 

l(i 

-( 2 )^ + 

„ xx 
(m-r 


kaby ' 


S 22 ~ _S 1 1 


(5.20) 


where we have defined 


and 


f = (a - Y f x )/ a > m = - Y f y / b 


2 2^/2 

v = (1 - I - m ) 


(5.21) 


(5.22) 


and where again we have 


2 2 
a =.l + f 

x 


b 2 = i + f 2 

y 


The wave numbers k4 and km are the projections of the propa- 
gation vector kj onto the tangent plane in the x and y directions 
respectively. The quantity v is a bit difficult to interpret, but it 
can be written in terms of the "tilt angles" of the surface, defined 


by tamL = f and tan 6 = f : 
7 T x y 



38 


~ cos ^ (0 - 40 - cos^0 sin^6 (5.23) 

At (locally) normal incidence, 4 J = 0 and 6 = 0 so that v = 1 
(its maximum value). For local incidence near grazing ( 4> = 0-ir/2) 
v may go to zero causing S to blow up. Since S must be small if 
it is to be a good approximation, large incidence angles must be 
avoided. 


The correction is seen to be identical to an asymp- 

totic form found by Fock (1945) for the current distribution on a 
convex paraboloid of revolution. Fock did not use the integral 
equation, but rather solved Maxwell's equations by separation of 
variables. Fock's asymptotic form holds for large distances away 
from the shadow boundary. 

Jackson (1971) has shown that the corrected current density 
j(°) + jO) corresponds with Rice's current density in the case of 
gentle curvature and moderate incidence angles. 

It is interesting to compare the criteria (inequalities (2.7)) 
with our correction. For the one - dimens ional case, f = f(x), alone, 
say, the condition that S be small gives 


xx 1 


<< 1 . 


2kv~ 


3 - 1 

Or, since | f ^ | /a = r^ , v = cos(0 - \\j) - cos 0 1 , this means that 


4-n-r cos 0 1 > > \ 
c 


(5.24) 


Note that S is purely imaginary in number. This means 
that the perturbed current density is 90° out of phase with the 

zeroth-order current. Exactly how this phase shift will 
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determine the scattered field will depend on the he ight- curvature cor- 
relation properties of the surface. 

Returning to the point made earlier, that there Is an ap- 
parent contradiction in the integral equation between the local nature 
of the high frequency approximation and the importance of far-source 
contributions, we examine the integral S (see the Appendix). The 
inductive contribution comes from the integrals 

oo -ikp.. 

Ind. = f J n (Lp^)e dp^ 9 n = 0, 2 

o 

and the radiative contribution comes from the integrals, 

oo -ikp^ * 

Rad-” | ikpj J (Lp^)e dp^ , n = 0,2 . 

o 

-2 - 1-2 

The ratio Ind/Rad = Rad = v for n = 0 and Zv + v for n = 2* 

±hus, inductive and radiative contributions are comparable: this 
despite an intuition that the inductive component might preponderate* 

3 

The rapid increase of the ratio with incidence angle (v oc sec 0) 
warns us that far-source contributions are becoming more important, 
and that the surface curvature must be increasingly mild if the local 
nature of the solution is not to be violated. 

To summarize our results, we write the corrected current 
density J. as 

j = j ( o) +2 (1) = (i + §) J (o) (5.25) 

where 1 is the identity matrix* 
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VI. Calculation of the Scattered Power 


We apply our results to the calculation of the average power 
backscattered from a random rough surface. In the case of perfect 
conductivity, the far-field Stratton-Chu integral'' reduces to: 


-ik e 
4ttR 


R X J/Rx(-r|J secure 

A S 

o 


E_ is the electric field vector, R- RR is the position vector of the 
far field (Fraunhofer zone) point; x is the position vector of the source 
points on the surface; is the illuminated area; and r\ is the < 

impedance of free space. 

For backscatter, the unit vector R is directed toward the 
source of incident radiation and so is given by 


R = (-a, 0, y) . 

Since in the far field the E -vector oscillates transversely to the propa- 
gation vector kR , we have 

E -R = 0 ; 



For example, see Silver (1949), p. 161. The Stratton-Chu integral 
is a particular vector form of the Helmholtz integral, eq, (2.1). 
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n J sec co = 
— s 


y 

j j + f j , 

x x y y ± 


and J and J are given by equation (5.25). (We are now assuming a 

x y 

unit incident electric vector.) Expressing the amplitude vector of 
the incident field A explicitly for H and V polarizations, but keeping 
the symbols S - j for the elements of S ^ equation (6.1) yields for the 
four polarization combinations : 


E HH = CffUy + a £ x )(l - Sj L ) - a f y S 2 j] e l2k ( ax “\ f ) dxdy- (6.2a) 

E HV = -C//[(y + af/s i2 - 2af y (y + af y )S n 

-(af y ) 2 S 21 ]e' l2k * aX '^dx dy (6.2b) 

E VH = CffS zl e- i2k(aX -y £ ) dx dy (6.2c) 

E VV = - C//[( y+af x )(l+S 11 ) +af y S 21 ] e " l2k ( ax -V f ) dxdy (6. 2d) 


The f irst H or V stands for a horizontally or vertically polarized 
incident wave; the second H or V stands for the horizontal or 
vertical component of the backscattered wave. In the above, we have 
used the fact that = ~S ^ ; we have let C stand for 

- ik exp(- ikR) / 27rR . 

2 

The scattered power is proportional to | E | . The usual way 

2 

to calculate |E| is to form the two-fold integral from |E| = E E r . 
The integrals (6.2) are of the form 
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E = CJ"/(K + P) e _l 2 k(ax 'V f ) dxdy (6.3) 

where K is the Kirchoff term, equal to y + a for co -polar ization 
(IfHjVV) and equal to zero for cros s -polar ization (HV, VH). P is 
the perturbation part. The forms (6.3) yield for the co-polarized EE' 1 ': 


EE* = CC*ffff(KK' + K’P^+KP' +P 1 P >:: )e" l2k ^ a(x ' x > ' v( f ' - f )] dx dy dx < dy t . (6.4a) 
and for the c ro s s -polarized EE * : 


EE*= CC*//JjP'P s V l2k ^ X '“ X * v( J -'“ f )] dxdydx' dy 1 


(6.4b) 


Jackson (1971) discarded the P l P '' terms in co-polarized re- 
turns on the grounds that their magnitude was of the same order as the 

2 Z 

error in the cross-product terms, L e. , of 0(\ /r^ )• This is a reason- 
able choice for small angles of incidence where the K terms are of 
0(1). Toward high angles of incidence, however, V + a becomes 
appreciably less than unity, and can have a magnitude O(f^) = OfX/r^). 
Thus, at large incidence angles the error in the cross-product terms 
may be significantly less than the magnitude of the P ! P* terms, justify- 
ing the retention of the P’P' 1 ' terms. Also, we must consider the fact 
that we do not know how the errors will transform through the integral 
operation. In any event, we should not want to lose possibly valid 
information, and moreover, we should like a guarantee of the power 
being a positive quantity by keeping a completed square. 

We consider a rough surface z = f(x, y) to be a realization of 
a stationary (homogeneous) random process. The return powers are 
then random variables. For an illuminated area very much larger 
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than the scale of roughness, one might expect the variability in return 
power to be small. However, this ignores the fact that scattered 
radiation from different portions of A q will have random relative phases, 
resulting with Rayleigh-type statistics of signal fading and reinforce- 
ment. To calculate the average power we proceed by taking an ensemble 
average or expectation of all possible surface realizations, denoted by 
corner brackets <(*•*)>. Since expectation and integral operations are 
commutative, the average power return in (6.4a) can be written as 

< | E | 2 ) = CC*ffff(K'K+ K'P*+ KP 1 + P'P*) • 

, • e i2k v( f '- f ) >e -i2ka(x'-x) dxdy dx , dy , f 

and similarly for (6.4b). 

An Immediate consequence of stationarity is that expectations 
of the type 

<K'P* + KP')e t2kY(f '' f) > 

can be expressed as 

$(|) + 

where 

$ = <K'P*e l2kY ^ f " f ^> 

and 

£ = x 1 - x . 

The expectations are computed on the assumption that f is a 

stationary Gaussian random process of zero mean, (f) = 0. Define 

the twelve dimensional random vector Y. whose first six elements 

are f. f . f , f , f , f and whose second six elements are 
x y xx xy yy 

f T , f^., etc. The mean of any derivative of a stationary process is 


zero; hence 
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< Y ) = 0 . 

Since the mean of the vector is zero the covariance matrix can 
be written as 

A(i) = (X Lj ) = ( YjYj) 

The multivariate Gaussian distribution with the covariance f\ has 
the probability density function, 

, , i , r .i t „-i 

p(i) = — t r exp i'2-^ a y 

(2tt) (detA) 2 

T • -1 

where y is the transpose of and A is the inverse of A* Define 
the characteristic function of JY , 

P(t) = (e 1 - — ) ; 

00 OO OO 

p (t) = / dy x / dy 2 * • • / dy 12 e - ' 1 p(^) . 

- 00 - OO - OO 

If p(^r) is the multivariate normal distribution (6.7), then has the 
form"' : 

<p(t) = exp j t_ T Ai 1 • 

Or, in terms of elements L , 

y (tj. fc 2* * • * > fc i2^ = exp {" 2 r X ij t L t j 1 

Now, the required expectations can be generated in a simple 
manner from the characteristic function by expressing the slope- 
dependent coefficients in S as polynomials in f^ and f ^ . We expand 
the (three) slope -coefficients in (5.20) in a Taylor series about 

*E. g. , Wilkg (1962), p. 168. 



( 6 . 6 ) 

(6.7) 

(6.8a) 

(6.8b) 


(6.9) 
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f x = f . = 0. (Expansion about the rms values \f\ ^ , \j\^ ^ might be 
more sensible, but it is a good bit more difficult.) We can truncate at 
first, second, or third order in slope. With the multipliers of the 
phasor exp { i2ky(f ' - f ) } expressed as polynomials in the Y-elements 
we can compute term by term the expectations of the forms: 

i2ky(Y 7 - Yi) 


<e 


) 


i2ky(Y -YJ 

(Ye > 

P 

i2ky(Y -Y ) 

<Y p Y q e > 


( 6 . 10 ) 


These averages are computed from the characteristic function in the 
manner outlined: Define the twelve dimensional vector t' v (star does 
not mean complex conjugate) all of whose elements are zero except for 


t^' and which have the values, 

t L * = -2ky 


t^ = +2ky 


( 6 . 11 ) 


The expectations (6.10) can then be written as 

/ Y \ 

(e- — > 


(Y e - ) 

N P 

(Y D Y a e - ' - ) 

p q 


From the definition of the characteristic function (6.8), we find 
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(e l ~ ' - > = <Mt*) 


p q 


-i a* 


at 


p 

t_= t'* 

_ .-2 

A 1 

at at | 


p q 


t = t" 


( 6 . 12 ) 


From (6.9) and the definition of t’ 1 ' (6.11) we find 


2 2 , 


■Mtf) = exp {-4k v (x n - x 17 )} 

tt = 2k vV - V * ( - 1 


d<j> 

at" 


d 2 <p 

at at 
p q 


t* 


[4k 2 7 2 (-X ql + X q7 )(-X pl - + X p7 >- 


(6.13) 


In the manner outlined, the required expectations can be 
calculated. The last step is to find the covariances . as a function 

y 

of the lag ^ = x 1 - x* All 7Z covariances can be expressed as partial 
derivatives of the covariance function, 

R(D = (ff). 

In accordance with our "smoothness 11 condition, R(|) possesses con- 
tinuous partial derivatives of all orders. 

If the illuminated area is large compared to the scales 
of roughness in the x and y directions ("correlation lengths"), then 
the scattering integrals of the form 



dxdydx'dy* , 


k_ 

2ttR 


//// C(e)e- i2kQ ( x '- x > 


are nearly equal to 



7 C(i) e _ L ^ ka ^ d| dr| . 
-Y 


C(£) stands for the expectation ( {• • • }e* 2 ^^ ^ ) and Y is a large 

distance in the rj direction. Because ot the behavior ot the exponential 
< p , Y can go to infinity only in the sense of the limit 

lim / / C (|) e " l2ka ^ d£ dp ; 

Y— oo -oo -Y 


and it is in this limiting sense that the infinite limits oi integration 
in the final formulas have been applied, v 

The power returns are usually given in terms of the normal- 
ized isotropic radar backscatter cross sections, <r° , defined by 


<r° 


4ttR 2 

A 


<|E| 2 > . 


o 

2 

< | E | ) is the quantity we have calculated, namely, the ratio of 

2 2 

(|E| ) at the receiver to |E| incident. Here, we do not consider 

realistic antenna gain patterns. The incident field is taken to be of 
constant amplitude over the area A Q .' r > ' 


* Cf. pp. 13, 14. 

-"Requiring the incident field to be of constant amplitude is only an 

artificial restriction. The cross-section is really determined by 

the statistical structure of a small patch, over which the antenna 

gain pattern is practically constant. The effect of gain pattern 

can be accounted for formally by introducing an effective area 

A in place of A . 
e L o 
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Following are the scattering integrals arrived at. The co- 
polarized returns are expressed as 


°HH 


= K + P + S 


c r 


VV 


= K - P + S 


(6.14) 


where K represents the "Kirchoff" part, P the M perturbat ion" part, 
and S represents the perturbation "squared". The P-term is calcu- 
lated to* second order in slope, and the S-term is calculated to zeroth 
order in slope. To first order in slope, the cross -polar ized returns 
are equal: 


cr - cr 

HV VH 


o 



(6.15) 


where D stands for "depolarized" return. As previously defined, 
B stands for 

B - -R(0) + R(|) ; 

we introduce the additional symbols: 

* etc * 


We find, then, that: 
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£°o oo , ? _ 72 7 

K = — jj— / d£ | dr> J [ -y - a {R|^ + 4k ^ R ^ )] cos ’.Ka| 

O - oo l 

- tkY^aR^ sin 2ka|j- 


,4k 2 y 2 B 


(16a) 


2 oo oo - 2 2 2 2 2 2 

P = — — J d£ / dti |l2(B|| - v B^) - 8k a B ||( R | ‘ Y R T) ^ 


2 -2^2 


2-2 


+ 4a V ' B || " 2q Y " Y 


+ 4a 2 Y _2 (R^ R ^|| - Y^ R ^££^1 cos 2ka t 

+ [4ka(R ^ B|| - +Y 2R e B T1 T 1 ) 


-1 


-k'ay “(R^^ - y“ r 1 ^] sin 2ka 1- e' 11 ' Y 




2 2 
4k v B 


(16b) 


s = ^/d5 /d„j v -‘<B 4e - v“B^) 


Z 0O oo 


f -2, 


V 2 "4 
k V ■ 


(r £US 


- 2 v R 


ItTR] 


4„ si ,, t 4k 2 Y 2 B 

+ v R )y.cos2ka|e 
Y •n'n'n-n' r 


(16c) 


°d = / dr > Tt 4B lr, + M 2r 64t)^ cos 2ka ^ 

o - oo L 

+ [l6kaB |Ti (R | B ( . Ti + R T i B U ) 

i ? n . ,'\ 4k 2 v 2 E 

+ ik- 1 Y - 2 a(B w R e ^+B^R WT| + R t| R mTI + R fc R feetltI )l«m2ka6|e 


(I6d) 
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The Kirchoff term is in the form used by Chia (1968). Beck- 
mann's form, eq. (2.22), differs because of the integration by parts 
and the neglect of edge terms. The K-term can be integrated by 
parts, yielding Beckmann's form upon neglecting "edge terms". For 
most practical applications, we have deep fade conditions, 2ky(r > > 1, 
and a large A q . One need only be careful of edge effects when 
numerically computing the integral for large angles approaching 
grazing incidence (this would entail having to carry the integration 
over a larger lag). 
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VII. Prop er ties of t he Solution 

Here we examine the general nature of the solution, give 
a sample calculation, and present for comparison some published sea- 
return data. 

7.1 Compariso n with the results of first-or d er Rayleigh-R ice theory 
It was shown in section 3.3 that for values of the spectral 
constant A typical of the sea, the small perturbation approximation be- 
comes very accurate toward large incidence angles. This fact demon- 
strates the validity of small-perturbation methods for sea-return from 
angles of incidence exceeding 0 = 45° or so. We should like to compare 
our solution with Rayleigh-Rice theory for the condition of small rough- 
ness amplitude (i. e. , for 2k^cr<< 1). If we linearize the integrands of 
the scattering integrals (6.16) and integrate using the definition of the 
spectrum (eq. (3. 12) ), we get scattering c ros s -sections of the form: 

0-° = 8irk 4 g(0) S(2ka, 0) 

where 

g HH ( 0 ) = (l-ta n Z 8) Z (7.1a) 

g vy (0) = (1 +tan 2 0) 2 (7.1b) 

and the depolarized return is zero. Comparing this with Rice's 
g(9), eqs. (4.4 ), viz. , 

4 2 2 

g HH (9) = cos 0 = (1 - sin 0) 

gyy(9) = (1 + sin 2 0) 2 

we see that (7.1) is in good agreement for 0 ^30° . The failure to 
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match Rice’s g(0) for large incidence angles indicates a general failure 

of our solution for representing backscatter from any surface having an 

* 

appreciable spectral density at the first-order Bragg condition (K = 2ka) 
for 0 ~ 40 . Clearly, the simple curvature correction is over -predicting 
the splitting of the copolarized returns from these angles. The failure 
of representation at large angles is not surprising, and could have been 
expected on the basis of the criterion (5.24). 

» 

1 . -4 

7 . 2 A sample calculation : Copolarized returns from an isotropic K - 

o 

spectral-law surface and a one-dimensional K~ J spectral-law surface 

We should like to apply our formulas to scattering from an iso- 

,-4 

tropic surface described by the simple K spectral law (eq. (3. 17)). 

-4 . 

It is understood that a K law surface is extremely rough (having 

first-order discontinuities) and poorly conforms with our requirements 

-4 

of smoothness. However, as the K law spectrum is a simple and 

realistic descriptor of the sea surface, we shall use it. 

-4 

Assume that the K law is vaiid for wave numbers in the re- 
gion of the Bragg condition K = 2ka. The detailed behavior of the 
spectrum for wave numbers much smaller than or much larger than the 
Bragg wave number little affects the scattering. As in section 3.3 
we take a low wave number cutoff corresponding to the spectral 

peak (cf. eq. (3. 18)). The covariance function is very well approxi- 
mated by eq. (3. 19)* Now, the second derivatives of the covariance 
function diverge logarithmically at the origin. As far as the problem 
of computing the scattering integrals is concerned, we can get around 



this mathematical difficulty. We can impose a high-frequency spectral 
cutoff, K = , say, corresponding to a viscous cutoff. Cutting off 

the spectrum abruptly this way creates a high frequency oscillation 
near the origin which is entirely artificial. Rather than doing this, 
it is simpler (and more sensible) just to "cut out" a small lag region 
near the origin, r ^r^, in which the second derivatives are to be 
taken as constant. Thus, we define for r <r 


2 9 R ' 1 8R 2 

VRs -2 

8r 


2 . 


where c r is the slope variance. In terms of a high frequency 
s 

2 

spectral cutoff, o- would be given by 
s 

2 A . K „ 

% = Aln ir • 


(7.2) 


(7.3) 


From the definition of the covariance function (eq. (3.19)), viz. , 

2 2 

R (r) = o- + A (-1 + y + In K Q r/2) 


(7.4) 


it is seen that r relates to the cutoff K as 
v v 


r = 2(K e y ) -i (7.5) 

We consider only the zeroth order slope terms in the P-integral. 

On transforming the Kirchoff integral and the first two terms of the 
P-Lntegral to polar coordinates (r,<£) and integrating over <p , we 
find ( using the above definitions): 
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<r° = K ± P 

K = 2k y~ J J Q (Zkar)e ^ r dr 
o 

00 2 2 
2 i 1 2 r 4W v R 

P - 2k A J a X n — J ( 2 ka r)e 7 rdr ( 7 . 6 ) 

V 

Figure 7 shows the copolarized returns given by the integrals 

.3 

(7.6) . The values of the constants used are: A - 4.05 x 10 (from the 
Pierson and Moskowxitz frequency spectrum with A = a/2); wind speed 
U= 7. 5 msec k = 2 cm" * (X-band); = 20 cm * (a nominal value of the 
viscous cutoff). In the isotropic case, the perturbation (P-term) is 
zero at vertical incidence, as we can see by putting a = 0 in (7.6). The 
failure of the solution at angles approaching 0 = 45° is seen in the rapid 
falling away of the horizontally polarized returns as K + P goes to 
zero. This failure occurs in the manner described by the linearized 
solution (eqs . (7.1)) which gives a = 0. 

A similar treatment can be given to the two-dimensional scatter- 
ing analog. In the two-dimensional case, the covariance function is 
given by eq. (3.21), and Bgg is defined to be 

f 0 > i ^ 

B tt = { A ln(£/S v ) , (7,7) 

*In Neumann and Pierson (1966), pp. 349-352. 
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-1 

where in terms of K , ^ = (K e^) . The first term (zer oth-or der 

v v v 

in slope) of the P-integral becomes 

00 ^ _ 

2 2 

P = 4kA I' ln(£/£ ) co s 2ka£ e 4k Y B d£ (7.8) 

e ' v 

V 

Fig. 8 shows the copolarized cros s - section , cr 0 = K _+ P, with P 

given by (7. 8) and K given by eq. (3. 7). The constants used were 

A=5 X]0 U = ]5m sec \ k = 2 cm * and K = 20 cm ^ . The values 

v 

of A, U, and k are the same values used for the Kirchoff integral shown 

by curve ’a 1 in Fig. 3. Again, the solution fails in the neighborhood of 

0 = 45°. The crqss-over of the HH and VV curves in the neighborhood 
° 

of 0 = 20 is an interesting result of the two-dimensional calculation. This 
cross-over is clearly evident in the 1.3 GHz backscatter data of McDonald 
(1 956) and the 8 . 9 GHz data of Daley et al. (1971), which data are repro- 
duced in Fig. 9. On the basis of our theory, it would appear that cross- 
over in copolarized returns in Fig. 9 is due to anisotropic wave condi- 
tions in the high-frequency portion of the wave spectrum. 

A calculation of the depolarized return and the higher-order terms 
in the P and S integrals has been avoided. The kind of approximations 

we have made to the covariance function and its second derivatives using 
-4 

the K spectrum are not easily extended to the higher-order derivatives. 

In dealing with third derivatives, one could take the third derivative of the 

covariance function as it stands with the r ^ singularity at the origin and 

simply let the singularity be removed through the integration. Fourth 

-2 

derivatives have a singularity proportional to the r which is not removed 
in the integration. 
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Fig. 9. Ocean backscatter data showing cross-over of copolarized returns. 
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A detailed analysis of how the scattering integrals might be 
“best applied” to the radar sea-return problem is outside the scope 
of this work. For X-band radar sea return, an accurate spectral 
representation of the high-frequency gravity and capillary waves in- 
cluding a viscous cutoff is one possibility. Then, all required deri- 
vatives of the covariance function will exist and the scattering integrals 
could then be computed. For lower -frequency microwaves, a detailed 
modelling of the covariance function near the origin is not necessary. 
But it remains a problem how to control most reasonably the behavior 
of the higher-order derivatives near the origin. We see that the 
practical problem of using the scattering integrals is intimately linked 
with the problem of error in this high-frequency approximation to the 
scattering problem. 
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VIII. Discussion and Conclusions to be Drawn 

A few years ago, Prof. R. K. Moore pointed out the need for 
extending Kirchoff techniques beyond the tangent-plane approximation.* 
With the curvature correction to the tangent-plane approximation, 
this has in a good measure been accomplished. 

The sample calculations illustrate the potential strength of 
this ’’corrected Kirchoff theory” for predicting radar sea-return from 
the near- vertical (specular regime) and from the transition region 
between specular and diffraction regimes. Because of the greater 
analysis needed to reasonably model the fourth-order mixed partial 
derivative in the depolarized scattering integral, no sample calcu- 
lation has been given. A calculation of the depolarized signal should 
prove most interesting, for it is not masked by the large specular 
component as the copolarized signals are. As the depolarized return 
depends essentially on the nonlinear component of the source distri- 
bution, this scattering solution is in apeculiarly good position to 
model depolarization near the vertical. Comparison with the second- 
order Rayleigh-Rice solution will provide a basis for testing the two 
solutions over the full range of incidence angles. 

On entering the pure diffraction regime the solution deterior- 
ates rapidly due to the over-emphasis on curvature effects and the 
increasing importance of non-local diffraction processes. It is interest- 
ing to note that the tangent-plane approximation alone provides more 

*In his paper "Scattering from Rough Surfaces" delivered to the Aug. 

1969 General Assembly of the URSI, Ottawa, Canada. 
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reasonable predictions, although it cannot account for the polari- 
zation dependence. The excessive splitting of the copolarized re- 
turns at larger angles can be controlled by artificial means, and this 
possibility has not been explored. For example, the wave spectrum 
could be filtered according to a smoothness criterion in order to re- 
move the high frequency waves which contribute to an intolerably 
large surface curvature. The derivatives of the covariance function 
corresponding to the "smoothed" surface could then be used in the 
scattering integrals. In this way, it would be possible to control the 
extent of splitting, and to establish a better correspondence with reality. 
This route is, however, not phys ically appealing. 

Rather than manipulate the solution by artificial means, one 
should admit the failure of the solution and recognize that at the 
larger incidence angles small-perturbation methods become applicable. 
If the Rayle igh-R ice solution shown in Fig. 6 is compared with the 
corrected K irchoff. solution, Fig. 8 , there is seen a fairly con- 

tinuous transition between the two solutions in the neighborhood of 
45 degrees. The idea of forming a '‘composite 1 ' solution suggests it- 
self. If cr^ 0 represents the high-frequency cross-section and 0 '^° 
represents the low-frequency cross section,' 1 ' then a contiguous 
solution for all incidence angles might be obtained from a weighted 
addition of the two cross sections: for example, 

* Here tr ° might be calculated to first-order in slope, and corrected 
for finite conductivity by the relationship, eq, (3.9 ). The small 
perturbation cross-section might be a first- or second-order 
. Rayle igh-R ice solution for finite conductivity, possibly including 
the artifice of mean-plane tilting (Valenzuela, 1968; Guinard et 
ah, 1971). 
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cr° = W(0)<r 1 ° + (1 - W(0))<r 2 ° . 

The weighting function W would be near unity for small angles and 
fall fairly rapidly to near zero at some critical angle in the neighbor- 
hood of 45 degrees. The complementary weighting of the two cross- 
sections guarantees that if both cross-sections were identical (say, 
were both perfect solutions), then we should have cr° = cr^° = o* 2 ° , 

From the standpoint of electromagne ti c theory, it appears that 
between the high-frequency and small- perturbation approximations, 
the radar sea-return problem is virtually solved. This is, however, 
within the framework of Gaussian surface statistics. Clearly, we are 
at a point of theoretical development where a more accurate description 
of the sea-surface in terms of its statistical properties is needed. It 
seems to make little sense to calculate the small corrections to the 
Kirchoff value near the vertical when the Kirchoff value itself is likely 
to be in significant error in its Gaussian form. While the height 
distribution of ocean-surface gravity waves is to a very good first ap- 
proximation Gaussian, the skewness of the distribution caused by the 
nonlinearity of the wave motion becomes increasingly important when 
considering slope and curvature distributions, and the joint distri- 
butions of these variables. Longuet-Higgins (1963) in a remarkable 
paper has shown how higher-order surface statistics may be de- 

rived from the nonlinear hydrodynamical equations of motion for 
gravity waves. Longuet-Higgins' development in terms of character- 
istic functions and "cumulants" ties in closely with the methods we 
have employed to calculate the statistical averages in the scattering 
integrals, and one is tempted to think of a rather effortless extension 
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of Kirchoff techniques to the non-Gauss ian solution. The problem 
with a transition to non - Gauss ian statistics is that the deviation from 
normality of gravity waves is quite different from the deviation of 
capillary waves. For example, gravity waves exhibit a positively 
skewed distribution function, while the skewness of capillary waves 
can be negative. For decimeter radars,a modelling of gravity waves 
alone is possible; but for centimeter radars, both gravity and capil- 
lary waves are important, and a simple non-Gaussian model may 
be near- impos s ible to develop. 

A current research effort by Prof. W. J. Pierson here at New 
York University is aimed at providing a more detailed description of 
the high-frequency wave structure in terms of the wave-height spectral 
density (Pierson et al. , 1971). Open-ocean measurements of the high- 
frequency portion of the wave spectrum (characterized by "chop", 
"ripples" and capillary waves) are non-existent, although some statis- 
tical information is available (e. g. , the estimates of slope variance 
by Cox and Munk (1954)). Hess et al. (1969) have shown that labora- 
tory measurements of wind-generated waves (conducted in "wind- 
wave tunnels") can provide useful information on the open-ocean wave 
statistics. Along with certain field observations (e. g. , the near-to- 
shore observations of Leykin andRosenberg (1970)), laboratory 
studies are yielding the increased knowledge of the high wave -number 
portion of the spectrum necessary for more accurate radar sea-return 
predictions. 
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Studies concerned with the effect of surface roughness on the 
natural microwave emission characteristics of the sea require the bi- 
static scattering cross-section . In this problem, one is concerned 
with an absolute power level, and a relative decibel representation of 
the scattered power pattern is inappropriate. The polarization effects 
in the near - specular direction will appear small on a logar ithmic plot 
compared to the polarization effects at large angles from the specular 
direction, although the polarization effects (diffraction effects) in the 
two regimes are of comparable absolute value. Thus, any effort to ac- 
count for pure diffraction processes in the bistatic cross-section must 
model these processes in. the near - specular direction, as well as in the 
pure diffraction regime at angles far from the specular direction. Since 
this corrected Kir choff theory is particularly suited to the near - specular 
regime, an extension of the scattering integrals to the bistatic case should 
be considered. 


For a formulation of the natural emission problem, see Stogryn (1967b). 
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APPENDIX 


Evaluation of an Integral 

The Fourier integral, 


£ =h J‘ J^ (1) M (2) d^dr, 


(Al) 


involves three types of integrals, viz. 
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On transforming to the elliptical coordinates (pj,cp ) defined by 

t ^ >1 

£ = — cos cp 
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these integrals become 
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where 


and 
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0 = tan 
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Integrating over the polar angle we get 
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where J and J-> are Bessel functions. Evaluation of the radiation 
integrals, i. e. , those integrals involving the factor ikp^ are a bit 
tricky, and are not to be found in standard tables of Fourier /Bessel 
transforms. If the substitution s = ik is made, the integrals be- 
come formally Laplace transforms which can be found in tables. 

We find that 
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n = 0 
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for 

n = 2 


(A5) 

(A6) 


1 (A7 ) 


(A8) 
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2 2 

where v = (1 - L /k ) 


2 


The ratio of radiative to inductive contributions is for n = 0: 
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and for n = 2: 
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Substituting (A8) into (A7) yields for the forms (A2): 
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and 
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